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and dropping all terms proportional to non-negative powers or the 
logarithm of this mass gap. As an alternative we write the regulated 
sum over zero point energies in terms of the heat kernel and apply 
standard heat kernel subtractions with modifications which take the 
boundaries into account. Finally we discuss to what extent these 
procedures are equivalent to the usual renormalization condition that 
tadpoles vanish. 
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1 Introduction 



In the 1970's and 1980's the problem of how to compute quantum corrections 
to the mass of 1+1 dimensional solitons received wide attention |^. These 
soliton models could be made supersymmetric (susy), and various methods 
to compute the quantum mass were developed (see P] and references in [Q, 
and 0). It turned out that if one repeated for the susy kink exactly the 
same steps as performed in for the bosonic case, the results for the quan- 
tum mass of the susy kink depended both on the regularization method and 
the choice of boundary conditions. In 1997, in view of renewed interest in 
extended objects due to dualities in string theory, the issue of why different 
methods and different boundary conditions for the quantum fields gave dif- 
ferent answers for the quantum mass M (and the central charge Z in the susy 
cases) was reopened [^. It was found that two of the most prominent regu- 
larization schemes used for these problems, mode number cut-off and energy 
(=momentum) cut-off, gave different answers; imposing different boundary 
conditions yielded different answers as well. As we shall review below, the 
reasons for these differences are by now well-understood, and both schemes 
yield the same correct answers provided they are applied correctly. 

In this paper we apply another prominent regularization scheme to this 
problem: zeta function regularization 0. Our aim is not to settle the ques- 
tion what the correct mass of the kink or susy kink is, but rather to show in a 
clear and simple example how to use zeta function regularization in problems 
with boundaries. We start from the regulated sum X] i^n~^* where s is a real 
positive number. To remove afterwards the pole 1/s we use two subtraction 
schemes often employed in non-renormalizable field theories (note however 
that here we are dealing with a renormalizable theory). The first scheme 
comes from Casimir energy calculations ^j. One just subtracts the singular- 
ities (including the constant term) in the limit of large mass (to be discussed 
below). The second scheme, widely used in quantum calculations on curved 
backgrounds, requires subtraction of the contributions from the first sev- 
eral heat kernel coefficients 0. The heat kernel formalism has been used in 
quantum theory since works by Fock [T^ and Schwinger [|ll|]. Given recent 
advances of the technique |p!2| , which take boundaries into account, this 
scheme becomes a rather universal tool for quantum calculations in external 
fields. In this paper we demonstrate that the two schemes are equivalent: 
large mass asymptotics is the same as expansion in the proper time of the 
heat kernel. Moreover, in two-dimensional renormalizable field theories they 
are equivalent to the requirement that tadpoles are absent, as we prove in 
Appendix B. 

Closer study of various regularization schemes for the kink and the susy 
^For general reviews of the quantization of solitons see 0. 
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kink in the period after 1997 revealed that the difference in answers was due 
first of all to the fact that by evaluating sums over zero point energies one 
included boundary energies for most of the boundary conditions considered 
10, ^. For particular boundary conditions in the kink sector (the so-called 
twisted periodic and twisted antiperiodic boundary conditions) the bound- 
ary energy vanished, and for these conditions mode regularization gave the 
correct results [|14|. The realization that for the susy kink boundary energy 



should be subtracted suggested a refinement of this scheme: one should av- 
erage over (quartets of) boundary conditions, so that upon averaging the 
boundary energies would cancel [|1^. For the energy cut-off scheme, on the 



other hand, it was realized that instead of an abrupt cut-off at some maximal 
energy one should use a smooth cut-off, or the limit of a smooth cut-off, and 
this yields an extra term which resolved the problem with the quantum mass 
both of the ordinary and of the susy kink . The mass of the susy kink has 



also been calculated by scattering theory methods g, which were applied 
to the bosonic case earlier in |]T^. (In |T^ even a renormalization prescrip- 



tion based on zeta function regularization was given. In this paper we shall 
justify this prescription.) Also some other schemes for evaluating the sums 
over zero-point energies were developed which remove boundary energies, in 
particular a method based on first evaluating the 8/ dm derivative of the sum 
(where m is the meson mass) and then integrating this result w.r.t. m, using 
the renormalization condition M'^^\m = 0) = [p!4| . 



For the reader who may wonder whether there really is a problem in de- 
termining the quantum correction to the kink mass by evaluation of \h^n-i 
we give an example which should remove any doubt, and which also plays 
a role below. One can impose susy boundary conditions such that for every 
nonvanishing bosonic solution cj^ there is a corresponding fermionic solution 
with the same eigenvalue, uj(^ = cj,^ (there are no bosonic or fermionic solu- 
tions with exactly zero energy for the susy boundary conditions (^0]) ). The 
one-loop correction to the susy kink mass is given by 

M« = (E l^-^n - E ^^^nO - (E l^-^'f' - E l^-^'J'') + AM . (1) 

Namely, one subtracts the contributions from the trivial sector, denoted by 
(0), from those of the topological sector, and one must add a mass counter- 
term AM , the finite part of which is fixed by requiring the absence of 
tadpoles in flat space (or far from the kink). Clearly, all bosonic sums cancel 
corresponding fermionic sums, but as AM is nonvanishing and divergent, 
even in the susy case, this procedure gives a divergent answer for M'^^\ 
One can explicitly evaluate the energy density near the boundary for susy 
boundary conditions, and finds then that it is also divergent 0, giving as 
the boundary contribution AM — M^^\ where M^^^ is the correct value of 



the mass of the susy kink |T3]. Subtracting the boundary energy from the 
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total energy then gives the correct result. Thus the global answer actually 
is correct for the entire system with susy boundary conditions, but it is 
insufficient for obtaining the mass of the isolated kink. 

By directly evaluating the local energy density one can obtain the correct 
value for the quantum mass by integrating this density around the kink. The 
first such method was developed in where a higher- derivative action was 
used which preserves susy and is canonical (in the sense of no higher time 
derivatives). This scheme gave a result for the quantum mass which agreed 
with the results of Schonfeld the 8/ dm method [Q, and the phase shift 
methods H. For the susy sine-Gordon model, it also agrees with the results 



obtained from the Yang-Baxter equation [|1J]. Another scheme to evaluate 
the local energy density employs the concept of local mode regularization 
[19|. This method also yields the correct result. 

After this summary of previous work done on the calculations of the 
quantum mass of the susy kink, we come to a new development: the ap- 
plication of zeta-function regularization to this problem. One converts the 
sum i ^f^ra"^* i^^*^ integral by using a function which has poles in k 
corresponding to each solution One can also convert each term o;^"^* in 
the sum into an integral by using the gamma function T{s). By writing the 
sum of these integrals as an integral over the sum, one encounters the heat 
kernel, about which an enormous literature exists. We shall give a simple 
self-contained account directly applied to the case of the kink and susy kink. 

In section 2 we discuss the bosonic kink with Dirichlet boundary condi- 
tions and use the large mass subtraction scheme. In section 3 we consider the 
susy kink and use the heat kernel subtraction scheme. In section 4 we sum- 
marize the methods. In Appendix A we extend the large mass subtraction 
scheme of section 2 to Robin boundary conditions. As a by-product we derive 
a simple relation between the contributions to the vacuum energy from the 
boundaries and from the bound states. In Appendix B we show that these 
subtraction schemes are equivalent to requiring absence of tadpoles. 



2 Large mass method for the bosonic kink 

The principle of zeta-function regularization of a sum \ ^ujn (we set = 1) 
is to replace this sum by 

E''' = Wi^l + ^^'r'' + Wi-l + ^'f' ^ (2) 

^ n ^ i 

where ujn = \Jk'^ + 'm^ are the zero point energies for the continuum spectrum 
of the quantum fields in a kink background (discretized by putting the kink 
system in a box with length L), and —nl + m^ = o;^ j are the squared energies 
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of the discrete spectrum, namely the bound states and the zero mode. For 
the kink on an infinite interval there is one bound state and one zero mode. 
The real number s is positive and large enough to make the sum convergent, 
and /i^ is a small mass temporarily introduced for the treatment of the zero 
mode (it will be used to move the pole due to the zero mode away from 
the starting point of a branch cut). Later we set fi to zero, and consider 
the asymptotics of E'^^^ for m ^ oo while keeping all scattering data Ki 
and kn fixed (these scattering data may themselves depend on m). Our 
renormalization prescription will be to subtract all non-negative powers of m 
in this limit (together with In m terms)^. We shall see that this indeed gives 
a finite result for 

It is clear from (|^) that the zero mode (k^ = m^) does not contribute to 
For later use it is convenient to keep the zero mode in (^ since then 
some cancelations can be seen more easily. However, we could drop the zero 
mode and still get the same answer. Strictly speaking, there is not even a 
mode exactly at zero for Dirichlet boundary conditions. Sometimes, as in 
application of the heat kernel technique (see sec. 3), both sums in (^ must 
be considered on an equal footing. One then has to replace {—nf — m?)^ by 
(— K^ — m^)2~*. This step is allowed in the positive spectrum only. Therefore, 
in the next section the zero modes will be excluded explicitly before applying 
the zeta function regularization. 

We first convert the continuum sum into an integral 

i?i°it = ^/^(^^ + M^ + /.^)i-^^ln0(fc,m) . (3) 

The function (f){k,m) has zeros at k = kn such that k"^ + = uj"^. The 
integration contour runs anti-clockwise and consists of one branch at k = 
^k + ie, a second branch at A; = ^k — ie , and a small segment — e < '^k < e 
along the imaginary axis. 

As the mass parameter m may also appear in Ki and kn, the function 
(j){k, m) depends also on m. This m should be kept fixed in the asymptotic 
limit. Therefore, we have replaced in u'^ = k"^ + by M^. We also 
replace m? by in the sum over bound states. The asymptotic limit we 

^Inside the logarithm the mass m must be divided by some scale parameter ji. Our 
subtraction procedure does not depend on the choice of this scale parameter. Indeed, if 
one chooses a different scale jl the logarithm changes as ln(m//i) = ln(m//i) — h\{fl/jl). 
The term on the l.h.s. of this equation and the first term on the r.h.s. should be subtracted 
since they are logarithmic in m. The second term on the r.h.s. should be subtracted as 
well since it is proportional to 

^The dependence of higher order correlation functions (e.g., two-point functions) on 
the mass differs, of course, from that of due to the presence of extra dimensional 
parameters (external momenta). The extension of the large mass subtraction scheme to 
these correlation functions would require some extra effort. 
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need is M ^ oo. The values of oOn depend on the boundary conditions, and 
the only place where boundary conditions enter in the calculations is in the 
choice of (j). 

We consider the super symmetric Lagrangian 

£ = -1 (9,$) (d^^) - - + U') ^ (4) 

with polynomial superpotential 

The kink solution reads $ = (m/-\/2A) tanh(mx/2). Here m is the meson 
mass, and A is the (dimensionful) coupling constant. 

In this section we consider the bosonic case, and impose Dirichlet bound- 
ary conditions^ rj{—L/2) = ri{L/2) = on the scalar field fluctuations rj. 
Then (j){k,m) is given by 

(t){k,m) = sm{kL + 5{k)) (6) 
where the phase shift 5{k) is given by [|T], 

ol/c) = — 2 arctan — — . (7) 

Along the upper part of the contour we can approximate sin(A;L + S{k)) 
by — ^ exp{—ikL — i6 (k)) because the term with exp{ikL + i6{k)) vanishes as 
L — > oo. Along the lower part of the contour we retain ^exp{ikL + i5{k)). 
We obtain then 

EfL^'' = ^/||(A:^ + M^ + /.^)^^^[(.A:L-lne-^(^)) + (.A:L + lne^^(^))]. 



(8) 

The contribution from the third part of the contour is independent of the 
phase shift, and hence it is dropped. Although = is formally a solution of 
the equation m) = 0, the corresponding mode is identically zero, and so 
the contribution from the pole at = must be dropped; this is a peculiar 
property of Dirichlet boundary conditions.^ The terms proportional to L are 

^Robin boundary conditions are considered in Appendix A. 

^There is a careful way to exclude the pole at fc = 0. In passing from the sum over n 
in Eq. (^) to the integral over k in Eq. ^ one has the freedom to multiply m) by a 
function f{k) without zeros and poles inside the contour whereby the contour is initially 
chosen to intersect the real axis at A: > 0. This function may be adjusted to give the 
product (j)(k,m)f{k) a finite nonzero value at fc = 0. As a result there is no obstacle to 
moving the contour across A; = 0. The additional contribution to the energy resulting from 
f{k) is independent of the background (in the case with 5{k) given by Eq. one may 
take f{k) = 1/fc) and may be dropped. 
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independent of the phase shift and equal to the contribution from the trivial 
vacuum, and as usual we subtract them. 

To estimate the limit oo oi the integral, it is useful to make a 

Wick rotation. We replace S{k) by —5{—k) in the second term because then 
the exponentials with 6{k) become equal after the Wick rotation. The first 
part of the contour is rotated to the positive imaginary axis {k = in with 
real positive k) 

- L term = \ J ^/""^'^'V ' - /^^'^^f- Ine'^^^^)] . (9) 



The second part of the contour is rotated downwards {k = —in, again with 
real positive k) 

^Wf _ L term = - 1 ^e-'-i\-'\>^ - - ^x^]^^' A[in e-^^(-)] . (10) 



For K, > M the sum of both terms contains the factor 

_ e*^(|-^) + ^-Ml-s) = _2i cos TTS . (11) 

For K < M the two integrals (almost) cancel each other so that only contri- 
butions of the poles in 5{k) at imaginary k need to be studied further. 

E^J^^'' - L terms = -- 7 — co^tis{k^ - - ^^)-2~^[-i^S{iK)] 

2 J TT OK 

+pole contributions. (12) 
We now use an alternative form of the phase shift[^ 

m = 



, k — imll , k — im 
In + In 



k + im/2 k + im 



(13) 



It is clear from this representation that the k integral has poles at k = ±m/2 
and K = ±m. We now see the reason for adding /i^ to in (^: it avoids 
that the starting point of the cut at A; = iM coincides with the pole at 
k = im (which is due to the zero mode). Each of the two contours contains 
a half-circle around these poles, but by combining these two k contours from 
< K < M, one obtains the two residues at the poles. These two residues 
precisely cancel the terms with bound states in the original sum. At this 

This is a particular case of the more general formula valid for all reflectionless 
potentials. 
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point we can put /i = 0. Hence we have arrived at an expression for the 
quantum mass as an integral over only M < /t < oo, and the bound state 
contributions no longer appear explicitly.]^ 

The expression for the quantum energy reduces to 



2 i TT ^ ^ dn 

M 



d 



, K — m/2 , K — m 
In '— + In 



K + m/2 



K + m 



(14) 



We rewrite the expression in square brackets as a sum of the limit of this 
expression for k — >■ oo and a remainder 



In 



K 



m/2 , K 
' fin- 



m 



K + m/2 



K + m 



m 2m, K — m/2 , k 

) + (ln /o +l^^- 

K K K + mj I K + m 



^f^). (15) 



With this decomposition inserted into the integral, the first term is explicitly 
evaluated for nonvanishing s, whereas in the second term we may set s = 0. 
The first, s-dependent, contribution is given by 



^(1) 



1 r°° dn 



2 JM TT 
COS ITS 



COSTCSiK 



.3m 



477 



(3m) M 



-2s 



r(3/2-g)r(g) 



(16) 



while the second, s-independent, contribution yields 



e: 



(0) 

(2) 



M 



7T 



-t: — 



M2 



m 



m?' /A 



+ 



2m 



3m 



(17) 



For s ^ the term J*? contains a pole 1/s proportional to as well as 



'(1) 

finite contributions proportional to M° and In M. According to our prescrip- 
tion this entire term must be discarded. The term e'^^^ vanishes in the limit 
M — >■ cx) and therefore must be kept in its entirety. Direct evaluation of this 
term aX M = m yields the quantum mass of the bosonic kinkQ 



1 



m 



4v^ 27r^ 



;i8) 



This is the correct result 0, |^, 0] first found in |T|] . 

^^In other approaches one finds the quantum mass expressed in terms of only the energies 
of the bound states; see, e.g., JTof . 

^^It is convenient to represent the numerator of the 2nd term in square brackets in ( pT[ ) 
as 2m = — m + 3m and then combine these two parts with the 1st and the 3rd terms 
respectively. Each of the two integrals resulting from this procedure is convergent and can 
be easily calculated. 
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Although the terms in the square brackets in ( |T^ ) are related to two 
bound states on the background of the bosonic kink, these two terms do 
not describe contributions of from the modes of the discrete spectrum to the 
vacuum energy. The expression (|l^) has been obtained by a contour rotation. 
During this rotation, the pole terms (which are genuine contributions of the 
discrete spectrum) were cancelled by a part of the continuous spectrum. Note 
that both terms in ([l^ ) are divergent while contributions of the bound states 
to E^^^ are always finite and do not even require a regularization (see eq. (^. 



3 Heat kernel method for the susy kink 

We now consider the susy kink. The classical kink solution satisfies di(f)K + 
U = 0, and the background solution = (px and ip = is invariant under 
rigid susy with parameter e_, as0 

6ip- = -(9i0e_ + dt(j)e+ - f/e_ . (19) 

Note also S(j) = —ie+ip- + ie-ip^. There are two sets of supersymmetric 
boundary conditions for the quantum fields which are invariant under the e_ 
transformations. The first set reads 

V\dM = 0, llj+\dM = 0, 

(ai-f/'(0))^_|aM = O . (20) 

The second set is 

{di + U\4>))r]\9M = 0, V-|oM = 0, 

(9i + f/'(0))7A+|aAf = . (21) 

For the spinor field it is sufficient to impose the Dirichlet boundary con- 
ditions for one component. The Robin boundary conditions for the other 
component are then determined by the Dirac equation, which is equivalent 
to the following equations on the spinor field components 

{di + U')ilj+-do^- = , (22) 
(-^1 + U') ^_ + 9o^+ = . (23) 

If we require that the ip^ component satisfy the Dirichlet boundary condition 
(the 2nd equation in (PO)), equation (pSf) then yields the Robin boundary 
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condition for ip- (2nd line in (|20|)). The same is also true for (pi]). We need 
these boundary conditions for the other spinor component as well because 
the heat kernel method uses second-order differential operators, so that we 
must exclude spurious solutions^- 

To discuss the cancelation of the bosonic and fermionic contributions to 
the vacuum energy, it is convenient to introduce the operators 

Q± = 9, ± U' (24) 

and build from them the second-order differential operators 

A_ = g+g_ , A+ = g_g+ . (25) 

The eigenfrequencies of the bosonic fluctuations and those of ip+ are defined 
by the operator A+, but the eigenfrequencies for ip_ are defined by A_. The 
operators A+ and A_ satisfy the following intertwining relations 

A_g+ = g+A+, g_A_ = A+g_ . (26) 

Moreover, the boundary conditions ( PP] ) and (^Tj) are compatible with these 
relations: if ip+ is an eigenfunction of A+ which vanishes at the boundary, 
then Q+ip+ is an eigenfunction of A_ which obeys the second condition in 
(pop. Conversely, if is an eigenfunction of A_ which satisfies the second 
condition in (pO]) , then Q^ip- is an eigenfunction of A+ which vanishes at 
the boundary. Therefore, the two operators involved are isospectral (up to 
the zero modes which do not contribute to the vacuum energy). This shows 
a complete compensation of the regulated sum of the zero point energies in 
the bosonic and fermionic sector. 

The property of isospectrality is crucial in our discussion. Nonsupersym- 
metric boundary conditions in principle also can be handled with heat kernel 
methods, but the treatment is more complicated. 

Each operator A± can be represented as A± = df — m? — V±, where m? 
is chosen in such a way that V± vanishes as ±00. The potentials then 

read V± = U'"^ ± U"U — m?. From now on we drop the subscript ± where 
this cannot lead to confusion. Let k"^ + m? be the eigenvalues of A. In zeta 
functional regularization the regularized vacuum energy has the form 

E.= ~^Y.{kl + M-y-' , (27) 



^"^The susy invariance of (20) under e_ transformations is easy to prove; in particular 
the last condition in (|2^) transforms into the full (j) field equation which becomes d/dt 
independent at the boundary due to the condition 77 = 0, and therefore this field equation 
itself is a boundary condition |Q. The first of the second set of boundary conditions 
should read £?i0 + U = Q. Then the first condition transforms into the second, the second 
into the third, and the third into {di + U')fj = 0. The last condition follows from the first 
by taking the time derivative. 
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where jl is an arbitrary parameter with the dimension of a mass which is 
needed to restore the correct dimension of Eq. Again, we have replaced the 
mass m in the spectrum by M to be able to consider later the limit M — cxd 
without deforming the potential V. The mass parameter m which enters the 
susy boundary conditions (|20| ) and (|2T| ) through the superpotential U also 



must be kept fixed in the limit M oo. The values of fc„ do not depend on 
M and are the same in the fermionic and bosonic sectors (up to the zero mode 
which must be excluded explicitly]^). Therefore, the bosonic contribution to 
the vacuum energy cancels the fermionic contribution for all values of M. 
Consequently, in the large M subtraction scheme both Eq and Eq" vanish 
for the susy kink. This gives the value zero for the total energy of the susy 
kink and the boundaries. As we shall see below, this is the correct result. 

To introduce the heat kernel subtraction scheme we first consider an indi- 
vidual contribution (bosonic or fermionic) to the regularized vacuum energy. 
We write 

i:(*^ + Mf-'^f . (28) 

where 

K{t) = e-"'' (29) 

n 

is the heat kernel for this problem. Its asymptotic expansion for t — > 0, 

K{t) ~ (30) 

n>0 

is equivalent to the expansion of Eq for M oo. The ultraviolet divergences 
are contained in the contributions with n = 0, |, 1 which correspond to the 
nonnegative powers of M. From this we define the divergent part of the 
vacuum energy as 



2 



n=0 



(aoF {s - 1) M2-2^ + ai r f s - M^-^^ + a^Y (s) M'^^ . (31) 

2r(s-i) I 2 V 27 J 

We define the renormalized vacuum energy as 

= Eq~ Ef^ (32) 

If one writes the zero mode contribution as square root of —rr? + M^, taking the 
hmit M — > cxD, the zero mode would contribute erroneously. 
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at s = and M = m. For n larger than unity, the coefficients 
multiphed by negative powers of M, so that the subtraction of the leading 
heat kernel asymptotics for small t is equivalent to the subtraction of the 
leading large mass asymptotics. In Appendix B it is demonstrated that for 
renormalizable theories these two procedures are equivalent to imposing the 
condition that tadpoles vanish. 

To obtain the kink mass, we must determine the contributions to the 
heat kernel from the boundary. This requires a detailed discussion of the 
heat kernel coefficientsl^. To write down these coefficients we need some new 
notations. Let us introduce the boundary operators B^'^ acting on a field (p 
(either bosonic or fermionic) 



(33) 



Here denotes partial derivative with respect to an inward pointing unit 
vector. The field ip satisfies the field equation {—d^ + V)ipn = ^"i^n- 

The asymptotic expansion of the heat kernel K{t) as t ^ for the bound- 
ary condition BifloM = (where B can be B^ or B^) reads 



E 



n=0,i,l, 



(34) 



with PI 



ao(i3«'^) 
a.{B^) 
ai{B^) 
ai(i3^) 



(An)- 
-I 

4 JdM 

(47r)- 
-(4vr) 



dx 



M 



dx 



ai{B 



- I Vdx + 2 

M JdM 



Vdx 



M 



1 

4 JdM 

Sdx 



dx 



(35) 



In the present (one-dimensional) case the boundary integral becomes a sum 
over the two boundary points. 

For the fermionic sector we sum the contributions from ip^ and ip^, and 
divide by a factor 2. It is easy to see that in the "total" heat kernel 



Jot 



„hos „ferm 



(36) 



the contribution of always cancels one half the contribution of rj. 

Let the kink be centered at the origin, but the box lie between x = Li 
and X = L2. In all super symmetric theories we have Oq"* = independently 

The expansion (|3^) has quite recently been used in numerical calculations of the kink 

mass 
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of the background and the boundary conditions. This is due to the equahty 
of the numbers of bosonic and fermionic degrees of freedom. The coefficient 
oi simply counts the number of Dirichlet minus Robin boundary conditions 
in each set (^) or (^TJ). Obviously, 

af = (37) 
for the boundary conditions (POf), and 



a'r' = - (38) 
for the boundary conditions (pi]). However, we must set a*i°* = for the 

2 

following reasonQ The operator A_|_ (resp. A_) has a zero mode which sat- 
isfies the Robin boundary condition Q+ip+ = (resp. Q-ip- = 0). Therefore, 
there is one zero mode in the fermionic sector for the boundary conditions 



(pop. (There are of course no zero modes in the bosonic sector because we 



use Dirichlet boundary conditions theref^ . Each bosonic (or fermionic) zero 
mode contributes +| (or — i) to the heat kernel coefficient. The zero mode 
structure, being a global characteristic of the problem, depends crucially on 
the boundary conditions (cf. |2^, |T5[ where the zero mode structure for 
other boundary conditions has been analysed). When the contributions from 
the boundaries to the vacuum energy are properly subtracted the mass of an 
isolated kink must not depend on the boundary conditions. In the zeta func- 
tion regularization, zero mode contributions must be explicitly subtracted 
from the vacuum energy. A similar mechanism works also for the boundary 
conditions (|2l|) . 

Let us now turn to calculation of ai. The volume terms in oi combine to 
give 



Jot 



1 f 

[vol] = -(47r)"^ / {-V+ + V-)dx 

2 JAI 

= -{Att)-^2 I {-2UU")dx 



1 



^(47r)-t(-2f/'(Li) + 2f/'(L2)) , (39) 



This follows also from the isospectrality arguments (see above). 
^^The zero mode ryo in the bosonic sector is a solution of the equation Q+i]q = 0. 
Explicitly, rjo = cdiipK, where c is a constant and (pK is the kink solution. On an infinite 
space the function -qo is indeed a normalizable zero mode. However, for finite L the 
Dirichlet boundary condition 77oU=±l/2 yields c — Q.A more detailed analysis shows that 
on a finite interval this mode is shifted and receives a positive energy. Note, that there 
was no zero mode contribution in the calculations of the previous section. 
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where we have used the Bogomolny equation 9i0 = —U{(f)). 

We next calculate contributions to oi from dM. We have to cast the 
boundary conditions into the standard form (pSl). Note that di = On at the 
left boundary, but di = —On at the right boundary. Consequently, we have 
to set S = —U' aX X = Li and S = U' at x = L2 for and the boundary 
conditions (pOl). For the boundary conditions (^) S assumes the opposite 
values, but now both rj and tp^ contribute. As a result, we obtain for (^Of) 
and (pID equal boundary contributions 

af [boundary] = ^{4n)-^2U'{Li) - 2U'{L2)) (40) 

which exactly cancel the volume term (^9|) giving 

af = . (41) 
We conclude that all singularities as well as the finite part of the contributions 



to the vacuum energy cancel. Thus, according to (^2]), not only £^0 but also 
the renormalized energy E'q™ of the kink plus boundaries vanishes. 

Qualitatively, it is clear what has happened. The vacuum energy associ- 
ated with the kink has been cancelled by the vacuum energy associated with 
the boundaries. Therefore, to obtain the correction to the mass of an isolated 
kink we have to subtract from the above vanishing result the Casimir energy 
of the boundaries as they are being moved to infinity. 

Of course, the kink solution cannot be smoothly deformed to the trivial 
one. We can, however, replace all background fields in the wave operator 
for the scalar field fluctuations and in the squared Dirac operator by their 
asymptotic values {V = Q), and replace U' in ( pi]) by +m at x = L/2 and 
—m at X = —L/2. In this way, we arrive at the problem of calculating the 
vacuum energy for free scalar and spinor fields with mass m. We consider the 
boundary conditions (^) because, as we shall later explain, the calculations 
with (pOD are more complicated. The boundary conditions for the free fields 
read 

{On - m)ri\oM = 0, ip-lsM = , 

{On - m)ij+\9M = , (42) 

where On = di on the left boundary x^ = —L/2, and On = — <9i at x^ = 
L/2. Obviously, the bosonic and fermionic contributions partially cancel 
each other. For ( ^2] ) the vacuum energy E^^'^ associated with the boundaries 
become 

= (43) 

where the superscripts [R] and [D] stand for Robin and Dirichlet boundary 
conditions respectively. 
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Let us start with the Robin sector. The eigenfunctions of the operator 

(pk = Ak sin(fcx) + Bk cos{kx) (44) 
satisfy the Robin boundary conditions in (^2]) if 

fkL\ (kL\ , , , 

cos I — 1 + m sin I — 1 =0 and -Bfc = (45) 

or if 

(kL\ ( kL\ 

— \ + m cos f — 1=0 and = . (46) 

It is easy to see that the two equations (|^) and (^) on the wave numbers 
k in are equivalent to a single equation, 

fR{k) = sin + 5i?(A;)^^ = with 6R{k) = arctan(/c/m) . (47) 

We can use the integral representation (|^) for the regularized vacuum energy 



E, 







where the contour again encircles the real positive semi-axis. 

Clearly, the same representation is also valid for E^^^ with Soik) = 0. 
The contour rotation proceeds in exactly the same way as before (although 
there are no bound states in this case). 

After rather elementary (but somewhat lengthy) calculations we obtain 
in the limit L ^ oo 

r.fBi. N 1 I /i^T(s)mi-2^ 1 , 

^'^'w- i^^FTTTiy"?-!^--- '''' 



Again we make the heat kernel expansion of the energy in the trivial 
sector (which is only located on the boundaries). Because the boundary 
conditions (which were susy in the kink sector) break susy in the trivial 
sector, both Eq and Eq^'' are nonvanishing in this case. From the equations 
( p5| ) we immediately obtain 

4^1=0, af = \, aP] = -^. (50) 



These coefficients define the divergent part of the vacuum energy (plf) . The 
boundary energy is the difference of the total result in the boundary energy 
in (^91) and the divergence in (^) 
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The one-loop correction M^^Ho the kink mass is the difference between the 
one-loop vacuum energy of the kink on a finite interval (which is zero) and 
the vacuum energy (51) associated with the boundaries. Thus we obtain 



M(^) = . (52) 

This is the correct result. 

Direct application of the same method to the other set of the boundary 
conditions ( pOD yields in addition to the eigenmodes (^) for ■?/'_ also the 
eigenmodes of the form sinh.{2kx / L) + Bf^ cosh.{2kx / L) which correspond 
to negative eigenvalues of —df. Their absolute values are large enough to 
make —dl + m? negative. In this case one cannot replace f/'(0) in the bound- 
ary conditions by its asymptotic values ±m before taking the limit L ^ oo. 
Keeping L-dependent boundary conditions considerably complicates the cal- 
culation. This case is considered from a little bit different point of view in 
Appendix A. 



4 Conclusions 

We have illustrated in a simple concrete physical model how one may use 
zeta function regularization and heat kernel methods. The problem was not 
totally trivial because there were boundaries present which played a role in 
the trivial and the topological sector. 

In the zeta function approach, one converts the sum of the modes into 
a complex contour integral. The renormalization condition which fixes the 
finite part of the zero point energy is that one must discard all terms pro- 
portional to non-negative powers of M or log M in the limit when the meson 
mass M tends to infinity while keeping the "scattering data" (namely ki and 
Ki) fixed. In our case these terms were proportional to and InM. The 
remainder gave the correct mass of the bosonic kink. 

In the approach which uses heat kernels, the sum over modes was con- 
verted into an integral over the proper time t, and the asymptotic expansion 
in t yielded the heat kernel coefficients a„. We first computed the total en- 
ergy in the susy kink sector, and defined its renormalization as subtraction 
of the contributions from ao, ai and ai. Both the total and renormalized 

2 

energy were found to vanish in the topological sector. This was to be ex- 
pected because we used susy boundary conditions in the kink sector, so that 
the sum over bosonic modes should exactly cancel the sum over fermionic 
modes. Then we repeated this calculation in the trivial sector, but with the 
same boundary conditions as for large x in the topological sector. These 
boundary conditions are not susy in the trivial sector. As a result we found 
that the total energy and the renormalized energy were nonvanishing. The 
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difference of the renormalized energy in the topological sector and in the triv- 
ial sector (which is thus minus the renormalized energy in the trivial sector 
in this case) is the kink mass. This difference yielded the correct result for 
the susy kink. A nontrivial technical aspect of this approach was the role of 
Robin boundary conditions and zero modes in the heat kernel expansion. 

The physical meaning of subtracting the contribution of the heat kernel 
coefficient is to subtract the energy of the trivial vacuum. However, the 
physical meaning of subtracting the contribution of the heat kernel coefficient 
ai is twofold: as expected it subtracts the tadpole contribution (the contri- 
bution from the counter term) and it adds the anomaly. It is well known that 
in scale invariant theories Oi contains the trace anomaly; here it contains the 
anomaly found in [|1^, ^. We shall discuss this elsewhere in more detail. 

Of course the problem confronting any global method using fixed bound- 
ary conditions in the kink and the trivial sectors, that there is a boundary 
energy for fermions in addition to the kink energy, also confronts heat kernel 
regularization. The method we use to overcome this problem, namely creat- 
ing a new problem in which only the boundary energy is present, presumably 
could be used in any regularization scheme. In particular, we have checked 
that it works for zeta function regularization with large-mass subtraction. 
We stress that the large mass subtraction scheme is fully equivalent to the 
heat kernel subtraction. It is possible, of course, to use the heat kernel sub- 
traction in the bosonic case and the large mass subtraction for the fermions. 
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Appendix A: Large mass subtraction scheme 
for Robin boundary conditions 

In this Appendix we demonstrate how the large M subtraction scheme works 
for generic Robin boundary conditions (Robin boundary conditions are a 
linear combination of Dirichlet and Neumann boundary conditions). All 
necessary information on one- dimensional scattering theory can be found in 
p3| . Consider a refiectionless potential V{x) locahzed around the point x = 
(and vanishing for x —>■ ±00). Before imposing any boundary condition we 
have two independent eigenf unctions (pi^2{k,x) of the operator A = df — 
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— V{x) for each k. Their asymptotics as x — ±00 are 

(pi{k,x) ^ cos{kx ± 6/2) , ip2{k,x) sm{kx ± 6/2) , (53) 

where 6{k) is the phase shift defined by the potential V. We put the system 
in the box —L/2 < x < L/2 and suppose that L is large enough so that near 
the boundaries the eigenfunctions can be represented by their asymptotics. 
Consider the Robin boundary conditions 

id. + ^i)¥^|(.=_^/2) = i-d. + ^2)¥^|(.=^/2) = (54) 

with arbitrary Si and 32- The general eigenmode of the operator A is a linear 
combination of ipi and ip2: ^{k, x) = Aipi{k, x) +B(p2{k, x). Substituting the 
asymptotics ( p3D of this eigenmode in the boundary condition ([5^ ) we arrive 
at the following equation which defines the spectrum of k: 

= (j){k) = sm{kL + 6 + 01 + 02) , (55) 

where 

ai,2 = — ai'ctan(/c/S'i,2) ■ (56) 

The vacuum energy is now defined by the equation (|^) where we have to 
substitute the function (j){k) given in (|55|) instead of (f){k, m). If we now repeat 
the same steps as in sec. 2 we arrive at the equation ( [T^ ) where 6 must be 
replaced hy 6 + a, a = ai + a2- The important property of the equation (|12[) 
is the linearity in 6. Therefore, the regularized vacuum energy for arbitrary 
auxiliary mass M can be represented as a sum of the contribution from the 
potential V (containing 6 under the integral) and the boundary contribution 
(containing a). To obtain the mass of the kink alone we have to subtract 
the whole boundary contribution, retaining the 5-term only. Hence we arrive 
at the old expression for the kink mass obtained for the Dirichlet boundary 
conditions. We conclude, that if the boundary contributions are subtracted 
the large M scheme gives equivalent results for the kink mass independently 
of the boundary conditions. 

We conclude this Appendix with a useful relation between the contribu- 
tions to the vacuum energy from the bound states and from the boundaries. 
For any reflectionless potential the phase shift can be represented as a sum 
over the bound states 



6{k) = 2 ^ In = 2 arctan ) , (57) 

k + XK^i, ^ \ k 



n 



where the bound state momenta. We can rewrite (|56|) in a similar 

form 



a = ai + a2 = arctan ( — j + arctan ) — vr . (58) 
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Because of the derivative dk in (|1^) the constant term — vr in (|58D does not 
contribute to the vacuum energy. By comparing (|57|) and (^8D we see that 
the vacuum energy associated with a boundary with Robin boundary con- 
dition^ is, roughly speaking, one half of the contribution of a bound state 
with the energy k = S. 

As an example, consider the supersymmetric boundary conditions (pO]). 
The bosonic sector is described by the potential and contains two bound 
states with the energies Ki = m and K2 = m/2. There is no boundary 
contribution to the vacuum energy. The same is true for ip^. The other 
spinor component ip- corresponds to the potential V- and has one bound 
state with K2 = m/2. For large L the Robin boundary conditions have 
Si = S2 = m. These two boundaries add the same contribution to the 
vacuum energy as would come from the bound state with ki = m. This 
restores the balance between the bosonic and fermionic contributions and 
gives zero total energy for the susy kink and the boundaries. A similar 
mechanism works for the boundary conditions (PTl), where 5*1 = 5*2 = — m 



and the boundary contribution is minus the contribution of the bound state 



Appendix B: Equivalence of the large mass and 
heat kernel subtraction schemes to the vanish- 
ing tadpole condition 

In this Appendix we demonstrate that in a renormalizable theory in two 
dimensions the heat kernel subtraction procedure is equivalent to imposing 
the renormalization condition of vanishing tadpoles in the trivial vacuum. 
As the former scheme is equivalent to the large M subtraction (see sec. 3), 
this will also prove that all three schemes are equivalent. For simplicity, we 
restrict ourselves to a bosonic theory. 

Consider first the case of a constant but arbitrary background field 
The second-order differential operator acting on the fluctuations has the form 
-dl + M{^f. For the action §) M^{^) = U'{^f + U{^)U"{^). The explicit 
form of will play no role. The complete heat kernel in this case is given 
by the term with n = in (^). As this term is proportional to the volume 
of the manifold, it is more convenient to consider the energy density Eo{x) 
which can be read off from the first term of the divergent part of the energy 

In the large M subtraction scheme the vacuum energy associated with the Dirichlet 
boundaries is zero. This is clear from the considerations of sec. 2. The reason is that the 
Dirichlet problem contains no dimensional parameter. Therefore, the boundary energy can 
only be proportional to M. In the large M scheme such contributions must be dropped. 



19 



Let us now expand Eq around s = 0: 



E,A.) ~ (; + 7. + 1 - In ^-^j - 9. Inr i^s - -) J - . (60) 

In a renormalizable theory ( 160D is accompanied by a counterterm 

A^;„(.) . _ (1 + + 1 _ _ a. ,„r (, - I) J Ml , (61) 

which repeats the pole structure of (|60D. Finite renormalization is encoded 
in an arbitrary parameter z/^ in (31). 



Vanishing of tadpoles means that the one- loop renormalized energy Eo{x) + 
AEq{x) has a minimum at the same value of $ = $o as the classical energy. 
Generically, d{M'^)/d^ ^ near $ = $o. We obtain: 

:iEoix) + AEoix)) = - In — + 1 = . (62) 



This condition defines z/^ in terms of M^($o)- 

Consider now a non-trivial vacuum $. The operator acting on the fluc- 
tuations is modified: —df + M($o)^ + V{x). It is essential that V{x) 
for X — > ±cxD. Consider ii^g'^ (|3lD. The first term containing Oq is now pro- 
portional to M($o)^; which is a constant. This term must be removed by 
renormalization of the "cosmological constant" . This is always done in such 
a way that the whole infinite volume contribution coming from oq is sub- 
tracted. The second term with ai is simply absent on manifolds without 
boundary. We are left with the term. It requires a counterterm 

- + ^ - (fl - ^ 4) J ^ ■ c'^) 

The fact the such a counterterm is indeed present is guaranteed by renormal- 
izability of the theory. As it is enough to do just one (mass) renormalization, 
the terms added to the pole - in the bracket in (|63D are exactly the same as 
in (|6ll) . Finally we have 

Et + AEo = -^(\n^ + l] . (64) 
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The right hand side of ( |64D vanishes due to the condition (0). We conclude 
that the renormalization condition that tadpoles vanish in the trivial vacuum 
is equivalent to dropping Eq" in the kink vacuum. 
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